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I INTRODUCTION 

In the last few years considerable experimental and theoretical work has been devoted 
to the study of quantum dots, which are atomic-like two-dimensional systems, confined 
within semiconductor heterostructures. * The parameters of such artificial atoms may dif- 
fer appreciately from their natural counterparts because of the interactions of the electrons 
with the crystal where they reside. In a quantum dot the natural atomic unit of length is 
a* = ea 2 /(m*e 2 ), where e is the dielectric constant and m* is the effective electron mass. 
Compared with the usual Bohr radius, ao = h 2 /(me 2 ), the length a* is typically large, 
e.g., a* ~ 185 ao in GaAs. The corresponding natural unit, with which we measure 
the magnetic field, B, is the field at which the magnetic length £b = he/ (S 1//2 c) equals a*, 
i.e., £?* = (ao/a*) 2 i?o, where B — e 3 m 2 c/?i 3 = 2.35 x 10 5 T is the value corresponding 
to free electrons. If a /a* is small, B* can be much smaller than B . Thus B* « 7 T in 
GaAs. This makes it possible to study in the laboratory effects which, for natural atoms, 
require the magnetic fields of white dwarfs or even neutron stars. 

The ground state properties of natural atoms in high magnetic fields have recently 
been analyzed rigorously in the asymptotic limit where the number of electrons and the 
nuclear charge are large [19-21]. For artificial atoms one may expect asymptotic analysis 
to be even more useful because the accuracy increases with the number of electrons, and 
a quantum dot can easily accommodate several hundred or even a thousand electrons. In 
the present paper we carry out such an analysis of the ground state of a quantum dot in 
a magnetic field. One of our conclusions is that the "self consistent model" introduced by 
McEuen et al [3,13] is a rigorous limit of quantum mechanics. This model has recently 
* The number of articles on this subject is by now quite large. See, e.g., [1], [2] for reviews, 
[3]- [7] for recent measurements of conductivity and capacity of quantum dots, and [8]- [18] for 
various theoretical aspects and further references. 
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been applied to explain interesting features of the addition spectra of large quantum dots 
in strong magnetic fields [15,6]. 

Before discussing our results for dots we summarize, for comparison, the main findings 
about atoms in [19-21] . The quantum mechanical ground state energy and electronic 
density of a natural atom or ion with electron number N and nuclear charge Z in a 
homogeneous magnetic field B can, in the limit N — > oo, Z — > oo with Z/N fixed, be 
described exactly by functionals of the density, or, in one case, of density matrices. There 
are five different functionals, depending on the way B varies with N as N — > oo. In 
each of the cases B <C iV 4 / 3 (with B measured in the natural unit Bq), B ~ iV 4 / 3 and 
iV 4 / 3 « B « N 3 the correct asymptotics is given by an appropriate functional of the 
semiclassical Thomas-Fermi type. For B ~ N 3 a novel type of functional, depending on 
density matrices is required, whereas the case B ^> N 3 is described by a density functional 
that can be minimized in closed form. 

A review of our results about quantum dots was given in [22]. Due to the reduced 
dimensionality of the electronic motion, there are only three different asymptotic theories 
for quantum dots instead of five for natural atoms. These three theories are given by 
simple functionals of the density and correspond respectively to the cases B <C N, B ~ N 
and B ^> N (B measured in units of £?*) as N — > oo with V/N fixed, where V is the 
attractive exterior potential that restricts the two-dimensional motion of the electrons. 
This potential, which plays the same role as the nuclear attraction in a natural atom, is 
generated in a quantum dot by exterior gates, and thus is adjustable to a certain extent. 
In the course of proving the asymptotic limits we shall also consider, in addition to the 
density functionals, a model of classical point charges in two dimension that gives a lower 
bound to the quantum-mechanical energy. 
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Some of the methods and results of the present paper contrast markedly with those of 
our earlier work [19-21]. From a mathematical point of view the most interesting feature 
of quantum dots compared to natural atoms is the somewhat peculiar electrostatics that 
appears because the interaction between the electrons is given by the three dimensional 
Coulomb potential although the motion is two dimensional. Also, the fact that the kinetic 
energy vanishes in the lowest Landau level requires additional mathematical effort in order 
to bound the kinetic energy from below by a functional of the density. We now describe 
in more detail the limit theorems to be proved in the sequel. A quantum dot with N 
electrons in a confining potential V and a homogeneous magnetic field B is modeled by 
the following Hamiltonian: 

H N = Y J H[ 3) + ^ £ \xi-xj\~\ (1.1) 

j = l l<i<j<N 

with Xi G R 2 and where Hi is the one-body Hamiltonian 

Hi = ^(iV-lA)\gJ^)s.B-(^)(^-\^)B + V( X ). (1.2) 
2m* V he J \2mc) \2mc ) \m* 2 J v ' v ' 

As before, e and m denote the charge and mass of a (free) electron, e is the dielectric 
constant, m* is the effective mass, and g* is the effective (7-factor. The magnetic vector 
potential is A(x) = \(-Bx 1 ,Bx 2 ) (with x = (x 1 ^ 2 ) E R 2 ), B = (0,0, B) and S is the 
vector of electron spin operators. The potential V(x) is supposed to be continuous and 
confining, which is to say that V(x) — * 00 as \x\ — > 00. It is not assumed to be circularly 
symmetric. The constant term in (1.2), —(he/(2mc))((m/m*) — \g*\/2)B, is included in 
order that the "kinetic energy" operator, H^ = Hi —V(x), has a spectrum starting at 
zero. The Hilbert space is that appropriate for fermions with spin, the antisymmetric 

N 

tensor product /\L 2 (R 2 ;C 2 ). 
1 
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We define an effective charge by e* = e/y/e and choose units such that Ti = m* = 
e* = 1. The unit of length is then the effective Bohr radius a* = 7i 2 /(m*e 2 ) and the 
unit of energy is E* = el/a* = e^m*/Ti 2 . Moreover, the unit £?* for the magnetic field 
is determined by heB*/(m*c) = E*, so £?* = e^m^c/(e 1 / 2 /i 3 ). The values for GaAs are 
a* = 9.8 nm, E* = 12 meV and B* = 6.7 T. 

The true quantum-mechanical ground state energy of H N is denoted by E Ci (N : B, V) 
and the true ground state electron density by p% B y( x )- The density functionals that 
describe the asymptotics of E® and p® are of three types. The first is a standard two- 
dimensional Thomas-Fermi energy functional 

£ TF [p;V] = (tt/2) J p(x) 2 dx + J V{x)p{x)dx + D(p, p) (1.3) 

with 

D( P , P )=\j j^Axly. (1.4) 

Here p is a nonnegative density on R 2 and all integrals are over R 2 unless otherwise stated. 
The second functional is a two-dimensional "magnetic Thomas Fermi functional" 

£ MTF [p;B,V] = J j B (p(x))dx + J V(x)p(x)dx + D(p,p) (1.5) 

where js is a piece-wise linear function that will be defined precisely in the next section. 
This functional is the two dimensional analogue of the three dimensional magnetic Thomas 
Fermi functional that was introduced in [23] and further studied in [24,25,21]. The present 
two dimensional version was first stated in [3]; these authors call it the self-consistent 
model (SC)*. 



* The repulsion term considered in [3] is slightly different from D(p 1 p), since it has cut-offs 
at long and short distances. It is still positive definite as a kernel and our methods can easily be 
adapted to prove Theorems 1.1 and 1.2 with such cut-off Coulomb kernels. 
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The last asymptotic functional will be called the "classical functional" , since the ki- 
netic energy term is absent and only classical interactions remain: 

£ C lP;V}= J V{x)p{x)dx + D{p,p). (1.6) 

The functionals (1.3) and (1.6) are in fact limiting cases of (1.5) for B — > and B — > oo 
respectively. As discussed in detail later, for each functional there is a unique density 
that minimizes it under the constraint f p = N. We denote these densities respectively 
by pj^ v (x), Pn T b,v( x ): Pn,v( x )i an d the corresponding minimal energies by E TF (N, V), 
E MTF (N, B, V) and E C (N, V). 

In order to relate E® to these other energies we take a high density limit. This is 
achieved by letting N tend to infinity (which is a reasonable thing to do physically, since N 
can be several hundred) and we let V tend to infinity. The latter statement means that we 
fix a potential v and set V = Nv. With this understanding of N, V — > oo our main results 
are summarized in the following two theorems. [In order to prove these theorems we need 
to assume that V is sufficiently regular. The technical requirement is that V belongs to 
the class C| '", (see Theorem 3.2 for the definition of Cj^'")]. 

1.1 THEOREM (Limit theorem for the energy). Let V = Nv with v a fixed function 



in C lo '". Then 



lim £ Q (iV, B, V) /E MTF (N, B, V) = 1 (1.7) 

N—too 



uniformly in B. Moreover, 



and 



lim E Q (N 7 B,V)/E TF (N,V) = 1 ifB/N^O (1.8) 



lim E Q (N,B,V)/E C (N,V) = 1 ifB/N^oo. (1.9) 
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1.2 THEOREM (Limit theorem for the density). Let V = Nv with v a fixed function 



inC\£. Then 



uniformly in B , and 



MTF /-, in \ 

Pn,B,V ~^ Pl,B/N,v V L - 1[J ) 



N 



JfPN,B,v^Pi?v ifB/N^O, (1.11) 

^Pn,b,v^P?,v ifB/N ^oo. (1.12) 

The convergence is in the weak L 1 sense* 

Let us add a few comments on these results. As discussed in the Section II, the energy 
£ MTF has the scaling property 

E MTF (N, B, V) = N 2 E MTF (1, B/N, V/N). (1.13) 

Thus (1.7) is equivalent to 

E Q (N, V, B) = N 2 E MTF (l, B/N, V/N) + o(N 2 ) (1.14) 

where the error term is uniformly bounded in B for V/N fixed. One expects the error to 
be (3(iV 3 / 2 ), which is the order of the exchange contribution to the Coulomb interaction, 
but our methods do not quite allow us to prove this. We do, however, show that for B/N 
larger than a critical value (depending on V/N) one has E MTF (N, £?, V) = E C (N, V) and 

£ Q (iV, V, B) > N 2 E c (l, V/N) - bN 3 / 2 (1.15) 

where the coefficient b depends only on V/N. 



* By definition, a sequence of functions f n converges to a function / in weak L 1 sense if 
/ fn9 ~~ > / f 9 for all bounded (measurable) functions g. 
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The condition that V/N is fixed as N — > oo guarantees that the diameter of the 
electronic density distribution stays bounded as N — > oo; thus the limit we are considering 
is really a high density limit rather than simply a large N limit. On the other hand, 
for a homogeneous potential V (e.g., quadratic, as is often assumed) one obtains also a 
nontrivial iV — > oo limit for V fixed, if the lengths are suitably scaled. In fact, this limit 
is given by the classical functional (1.6). Intuitively this is easy to understand, for if an 
increase in N is not compensated by an increase in V the charge density spreads out and 
the kinetic energy terms in (1.5) and (1.3) become negligible compared with the other 
terms. (The result again requires V to be in C lo '".) 

1.3 THEOREM (Energy limit with a homogeneous potential). Assume that v is ho- 
mogeneous of degree s > 1, i.e., 

v(Xx) = X s v(x). 

Then 

lim E Q (N, B, Kv) /E MTF (N, B, Kv) = 1 (1.16) 

AT— >oo 

uniformly in B and in K as long as K/N is bounded above. Moreover, if K/N — > as 
N — > oo, then 

lim E Q (N,B,Kv)/E c (N,B) = 1 (1.17) 

uniformly in B. 

One can also prove a limit theorem for the density in the case of homogeneous poten- 
tials. Since the formulation of such a theorem becomes somewhat complicated we refrain 
from doing this, but refer to Eqs. (2.14)-(2.16) below for the scaling of the MTF functional 
with k = K/N and to (3.24) for the weak coupling limit of the MTF density. 

The proof of the limit theorems involves the following steps. In sections II and III 
we discuss the basic properties of the functionals (1.3)-(1.6). In Section IV we consider 
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the energy of a system of classical point charged particles in R 2 in the exterior potential 
V as a function of the positions of the charges. This energy has a minimum, denoted 
by E P (N, V) (with "P" denoting "particle"). A significant remark is that the charge 
configuration, for which the minimum is obtained, is confined within a radius independent 
of the total charge N for fixed V/N. This "finite radius lemma", which also holds for the 
charge densities minimizing the functionals (1.3)-(1.6), is proved in an appendix. Using 
this and an electrostatics lemma of Lieb and Yau [26] we derive the bounds 

E C (N, V) - aN 3/2 > E P (N, V) > E C (N, V) - bN s/2 (1.18) 

where a and b depend only on V/N. These bounds are of independent interest apart 
from their role in the proof of the limit theorems where, in fact, only the latter inequality 
is needed. Upper and lower bounds to the quantum mechanical energy E Cl (N, B, V) in 
terms of E MTF (N 1 B,V) with controlled errors are derived in the final Section V. The 
upper bound is a straightforward variational calculation using magnetic coherent states 
in the same way as in [21]. For the lower bound one treats the cases of large B and 
small B separately. The estimate for large B is obtained by first noting that obviously 
E Cl {N,B,V) > E P (N,V), because the kinetic energy is nonnegative, and then using 
(1.18). For small B two auxiliary results are required: A generalization of the magnetic 
Lieb-Thirring equality considered in [21], and an estimate of the correlation energy. Once 
these have been established the proof of Theorem 1.1 is completed by coherent state 
analysis. The limit theorem for the density follows easily from the limit theorem for the 
energy by perturbing V with bounded functions. 

II. THE MTF THEORY: ITS DEFINITION AND PROPERTIES 

By employing the natural units defined in the introduction, the "kinetic energy" op- 
erator can be written 
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i^kin = \ (*V - A) 2 + 7 S • B - \ (1 - | 7 |) B (2.1) 
with 7 = ^m* /(2m). The spectrum of -f/kin is 

= (n + 7 (T+ (2.2) 

with n = 0, 1, . . ., a = ±1/2. We write the energy levels (2.2) in strictly increasing order 
as e u (B), v = 0,1, — The degeneracy of each level per unit area is d u (B) = B/(2tt), 
except if, by coincidence, 7 happens to be an integer; in that case d v {B) = B/(2tc) for 
v = 0, . . . , |7| — 1, while d u (B) is twice as large for the higher levels. It is worth recalling that 
if V(x) = K\x\ 2 , the spectrum of the one-body Hamiltonian Hi in (1.2) is solvable. The 
spectrum of Hi was determined by Fock [27] in 1928, two years before Landau's paper on 
the spectrum of (zV — A) 2 . For the Hamiltonian without spin, namely |(zV — A) 2 + K\x\ 2 , 
the spectrum is given by 

E = l(m - n 2 )B + l(m + n 2 + 1)[4K + B 2 } 1 / 2 

with toi, 77-2 = 0, 1,2,.. .. It is remarkable that this simple spectrum gives a qualitatively 
good fit to some of the data [8] . 

For a gas of noninteracting fermions with the energy spectrum (2.2) the energy density 
2b as a function of the particle density p is given by 
Jb(0) = 

MP) = £u(B) if D U (B) < p< D u+ i(B), z/ = 0,l... 
where j' B = djs / dp and 

D V {B)= J2 d -'( B )- 
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More explicitly, 



Jb(p) = £ e v {B)d v {B) + (p - D u _(B))e„_ +1 (B) 



1S = 



where z/ max = z/ max (p, B) is defined by 



(2.3) 



D^(B)<p<D„ max+1 (B). 

Thus js is a convex, piece- wise linear function with js(P) = for < p < di(B). It has 
the scaling property 

jB (p)=B 2 3l (p/B). (2.4) 
As B — > 0, becomes a quadratic function of the density: 



7T 



Moreover, 



for all p and £? (see Fig. 1). 



\imj B (p)=j (p) = -P 2 



3b(p)<3o(p) 




B/tt 2B/tt 3B/tt 4B/tt 

Fig. 1 The "kinetic" energy densities Jb{p) and jo(p) 
in the special case where 7 = 
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Given an exterior potential V the MTF functional is defined by (1.5). We assume 
that V is continuous* and tends to oo as \x\ — > oo. In particular, V is bounded below and, 
by adding a constant if necessary, we may assume that V(x) > everywhere. Because of 
(2.6) the functional (1.5) is defined for all nonnegative functions p such that J pV < oo, 
/ p 2 < oo and D{p, p) < oo. Since V > the functional is nonnegative. If N is some 
positive number we denote 

C N = {p\p>0, JpV < oo, fp 2 < oo, D(p, p) < oo, fp = N} 

and define the MTF energy by 

E MTF (N,B,V)= inf S[p;B,V]. (2.7) 

p6Cjv 

Because of (2.4) the energy satisfies the scaling relation 

E MTF (N, B, V) = N 2 E MTF (1, B/N, V/N). (2.8) 

In the limit B — > the kinetic energy density (2.3) converges to jo(p) = (ir/2)p 2 and 
(1.5) converges to the energy functional (1.3) of two-dimensional TF theory at B = 0. 
It is easy to see that also lim B ^ E MTF (N, B, V) = E TF (N, V), where E TF is defined in 
the same way as E MTF with (1.3) replacing (1.5). We can thus consider the TF theory 
as a special case of MTF theory. In the opposite limit, B — > oo, the kinetic energy 
term vanishes altogether and one obtains a classical electrostatic model (1.6) that we shall 
study in Section III. Note also that since j B < j for all B it follows that E MTF (N, B, V) < 
E TF (N,V) for all B. In particular £' MTF (1, j3, v) is uniformly bounded in the parameter 
13 = B/N for fixed v = V/N. 



V measurable and locally bounded would suffice 
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For fixed B and V, E MTF (N,B,V) is a convex, continuously differentiable function 
of iV and, since V > 0, it is monotonically increasing. By the methods of [28], [29], [21] 
(see also [30]) it is straightforward to prove the existence and uniqueness of a minimizer 
for the variational problem (2.7). 

2.1 THEOREM (Minimizer) . There is a unique density P^bV e ^ N suc ^ ^hat 

E MTF (N,B,V)=£ MTF (p™ TF v ). 

Note that the existence of a minimizing density with J p = N is guaranteed for all N 
because V(x) — > oo as |x| — > oo. This condition on V also implies that Pn,b,v vanishes 
outside a ball of finite radius, cf. Lemma Al in the Appendix. The scaling relation for the 
minimizing density is 

f%g,v(x) = Npffi NtV/N (x). (2.9) 

Theorem 2.1 includes the TF theory as a special case. In the same way as in Prop. 4.14 
in [21] one shows that P^ T J V — > Pjfv weakly in L 1 as B — > 0. 

The shape of the electronic density in the case of a quadratic potential V{x) = K\x\ 2 
and 7 = is shown in Fig. 2 for different values of B. The pictures were prepared by 
Kristinn Johnsen. At the highest value of B (8T), the density is everywhere below do(B) 
and given by the minimizer (3.15) of the classical functional (1.6). At B =7T, all the 
electrons are still in the lowest Landau level, but that level is full around the middle 
of the dot where the density is anchored at do(B). As the field gets weaker it becomes 
energetically favorable for electrons at the boundary of the dot, where the potential is 
high, to move into the next Landau level close to the minimum of the potential. A dome- 
shaped region then arises above the plateau at p = d (B) = D (B), but eventually the 
density hits the next plateau at p = D\(B). This gradual filling of Levels continues as 
the field strength goes down. At B = 2T three Landau levels are full and electrons in the 

13 



QDOTS EHLJPSJY 3 Fob /94 

central dome are beginning to occupy the fourth level. Finally, at B = 0, we have the 
usual Thomas-Fermi model, which may be regarded as a limiting case with infinitely many 
Landau levels occupied. 



zo 




8T 




7T 




20 



B 



OT 




Fig. 2 Quantum dots at various magnetic field strengths. The potential is V(x) = 
(l/2)m*a; 2 |x| 2 , with m* = 0.67m, %u = 3.37 meV, and N = 50. The coordinate axes 
are displayed in units of 10~ 8 m and the density p in the units 10 _14 m -2 . 



In order to state the variational equation for the minimization problem it is convenient 
to define the derivative j' B = dis/dp of the kinetic energy density everywhere, including 
points of discontinuity, as a set valued function (cf. [30]), namely 



{*u(B)} 



for D V (B) <P< D V+1 (B) , v = 0, 1, 



[e v (B),e v+1 (B)] for p = D V+1 {B), i/ = 0,l, 



(2.10) 
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With this notation the Thomas- Fermi equation for the functional (1.5) may be written as 
follows 

2.2 THEOREM (Thomas-Fermi equation). There is a nonnegative number fx = 
(jl(N, B, V) such that the minimizer p = PnzjV sa tisfies 

.-^-.•W-^fW' «%>J . (2.n, 

The quantity fx appearing in the TF equation is the physical chemical potential, i.e., 

p = dE(N,B,V)/dN. (2.12) 

Since E is convex as a function of N, fx is monotonically increasing with N for fixed B 
and V. It satisfies 

fi(N, B, V) = Nfji(l, B/N, V/N). (2.13) 

The derivation of the TF equation is analogous to that in [29]. It is also true that 
if (p, fx) is any solution pair for (2.11), then p is the minimizer of S for some N and 
fx = fx(N, B,V). The proof of this is a bit trickier than in the standard case [29], because 
2b is not continuously differentiable. It has been carried out by Lieb and Loss [30]. 

Finally we discuss the relationship between the MTF theory and the classical theory 
defined by the functional (1.6). We of course have that 



S MTF [p;B,V] = J 2 B (P(x))dx + S c [p;V]. 



From the definition of 2b one expects that the kinetic energy term above can be neglected 
for large B and hence that lims^oo E MTF = E c . The rigorous proof of this fact relies on 
a careful study of the classical problem. This analysis is far from trivial and is postponed 
to the next section. 
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There is another case where the MTF energy can be related to the classical energy. 
Namely, for a homogeneous exterior potential, i.e., 

V(Xx) = X s V(x) 

for all A > with some s > 0. We consider the potentials kV(x) with k > and are 
interested in the dependence of the MTF energy and density on the coupling constant k. 
Writing 

p( x ) = fe 2 /(«+i) /3 (j fe i/(«+i) a; ) (2.14) 

we have the scaling 



£ MTF { P ;B,kV] = k 2 '^) J Jb (p)+£ c [p;kV] 



= k 1/(s+1) 



^ 1/(s+1) Jj b (p)+S c [p;V] 



(2.15) 



where 



b = Bk- 2 ^ s+1 l (2.16) 



Changing k is thus equivalent to changing the kinetic energy by a multiplicative factor and 
rescaling the magnetic field, keeping the potential fixed. We shall show in the next section 
that for k small E is a good approximation to E MTF . 

III. THE CLASSICAL CONTINUOUS MODEL: A LIMIT OF MTF THEORY 

For densities p small enough (p(x) < d\{B) for all x) the kinetic energy js(P) vanishes. 
It is therefore natural to consider the resulting classical energy functional defined by 
(1.6), i.e., 

£ C [P-,V] = J p{x)V{x)dx+\ J J p(x)\x-y\- 1 p(y)dxdy. (3.1) 
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The corresponding classical energy is 

E C (N, V) = inf ^£ C [P; V] : P > 0, J p = iV j . (3.2) 

In this section we analyze this functional and prove that it is, indeed, the large B limit of 
MTF theory. 

As before we assume that the confining potential V is positive and that V(x) — > oo 
as |x| — > oo. Moreover, we shall here assume that V is continuous (in fact, we shall make 
an even more stringent regularity assumption in Theorem 3.2 below). 

We begin by showing the existence of a minimizer for (3.2). For general continuous 
V (without further assumptions) we must take into account the possibility that the min- 
imizing p may be a measure. In (3.2) we therefore minimize over all positive measures p 
with f p = N. It follows from the finite radius lemma given in the appendix that 



E C (N, V) = inf |^ C [P, V] : support p C {x : \x\ < R v }, p > 0, J p = Ar| . (3. 



3) 



Here R v depends only onu = V/N . Later on we shall show that the minimizer is, indeed, 
a function, and hence that (3.2) does give us the large B limit of MTF theory for suitable 
V. 

3.1. PROPOSITION (Existence and uniqueness of a minimizing measure). Let V 
be continuous. Then there is a unique positive measure p^ v with f P% y = N such that 
E c (N,V)=£ c lP% y ;V]. 

Proof: [Note: we write measures as p(x)dx, even if they are not absolutely continuous 
with respect to Lebesgue measure.] By (3.3) we can choose a sequence of positive measures, 
Pi, P2, ■ ■ ., supported in {x : \x\ < R v } with J p n = N, such that lim n ^ 00 £ c [p n ; V] = 
E C (N, V). The bounded measures are the dual of the continuous functions, and so, by the 
Banach-Alaoglu Theorem, we may assume (by possibly passing to a subsequence) that p n 
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converges weakly to a positive measure p still supported in {x : \x\ < R v }. In particular it 
follows that f p = N and f p n V — > f pV. Moreover, the product measure p n x p n — > p x p 
weakly. Hence 

y J P{x)p{y)\x-y\~ 1 dxdy=\uaJ J p(x)p(y)(\x - y\ + S^dxdy 

= lim lim / / p„(x)p„(y)(|a; - + ^"Maxh/ 

5^0 n^oo y y 

<liminf / / P n (x)p n (y) \x - yl^dxdy. (3.4) 

The first equality follows by the Lebesgue's monotone convergence Theorem. The last 
inequality is an immediate consequence of the pointwise bound (\x — y| + 5) _1 < \x — y\~ x - 
We conclude from (3.4) that 

E C (N, V) < S c [p- V] < liminf £ c [p n ; V] = E C (N, V) (3.5) 

and hence that p is a minimizer. 

The uniqueness of p follows from strict convexity of D(p, p). Q.E.D. 

The next theorem gives conditions, which are perfectly adequate for the physical 
applications, under which the minimizer is a function and not just a measure. Moreover, 
that function has certain nice integrability properties. 

3.2 THEOREM (The minimizer is a function). Assume that the potential V is in the 
class C lo '" for some < a < 1. (I.e., V is once continuously differ entiable and for each 
R> its derivative satisfies 

\VV(x)-VV(y)\<c R \x-y\ a , (3.6) 

inside the ball of radius R centered at the origin for some constant cr > 0.) Then the 
minimizing measure P% v of Prop. 3.1 is a function. It has the properties (with p being 
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the Fourier transform of p) 

j \p C N ,v(p)\ 2 \p\ r dp <C U -l<r<a (3.7) 
P% v * |a^| _1 is continuous (3-8) 
J p c Ny {xY&x < C 2 , l< q< ^— (3.9) 

where C\ and C 2 are constants (implicitly computed below) that depend only on the con- 
stants cr, q, r, a and on N. 

Proof: We write P% v = P- We know that J p = N and that p has compact support. 
From the former fact we conclude that p is well defined, continuous and bounded. 

Let g a be the function with Fourier transform 

UP) = J 9a (x)e^dx=^ (3.10) 

Then g a is continuous, / g a = 1 and / yg a (y)dy = 0. 
Let p a be the convolution p * g a , so that J p a = N. 

Since p is a minimizer, £ [p; V] < £ [p a ] V]. Explicitly this inequality is 

J(Pa-p)V + D( Pa , Pa ) - D(p, p) > 0, (3.11) 

Since f Vp a = f(V * g a )p we can rewrite the first term as 

J J[V(x-y) -V(x)]p(x)g a (y)dxdy. 

By integrating W along the line from x — y to x, and using (3.6), we have 

\V(x -y)- V(x) + y ■ VV(x)\ < c\y\ a+1 . 

(Note that by the finite radius Lemma Al all integrals are restricted to a finite ball.) Using 
that j yg a (y) = we can estimate the first term in (3.11) as follows: 



J (Pa -P)V<C j |y| a+1 |<7 a (y)|dj/ = (const. )a" 



a-l 
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The last two terms in (3.11) are 

(const.) J\p(p)\ 2 \p\- 1 [\g a (p)\ 2 -l]dp= (const.) J ^{p^W^p 



\p\>a 



(Recall that \x\~ x = ( const. 1 in two-dimensions.) The inequality (3.11) thus implies 

J |P(p)| 2 bl~ ld P < (const. )a" a_1 . (3.12) 

Using (3.12) and p(p) < J p = N we can now prove (3.7) as follows. 

OO 

J\p(p)\ 2 \p\ r dp = J \p(p)\ 2 \p\ r dp + J2 J \Kp)\ 2 \p\ r dp 

\p\<l n=0 2™<|p|<2™+ 1 

/oo „ 
|prdp + (const.)^2( r+1 )( n+1 ) / IP^^br'dp 



IpI<i 



n=0 



2™<|p|<2™+ 1 



< (const.)A^ 2 + (const.) ^2 (r+1 )( n+1 )- n ( Q+1 ) < oo, 

if r < a. 

To prove that p * is continuous is now easy. We simply prove that its Fourier 

transform is integrable. The Fourier transform of p * \x\ _1 is (const.)p(p)|p| _1 and 

/ 

J |P(p)lbl _1 dp< J |?(p)||p|- 1 dp+ 



VpI>i 





1/2 


( t V /2 


\p(p)\ 2 \p\ r dp 




/ \p\ r 2 dp < oo 


J 




\p\>i ) 



Finally, we prove (3.9). For 1 < q < 2 there is no problem because we know that 
J P 2 = J P 2 - Hence p is a square integrable function and, since J p = N, we conclude by 
Holder's inequality that (3.9) holds for 1 < q < 2. For q > 2 we will prove that 



/ 



pi* < oo with 



< t = 



< 2. 



a+ 2 ' q-1 

This will prove (3.9) by the Hausdorff- Young inequality, which states that (/ H*) 1 ^* > 
{jP q ) 1/q when 1 <t < 2. 
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We write = {|p(p)|*(l + |p|) m }{(l + m } and then use Holder's inequality 

with a -1 + 6 _1 = 1 to conclude that 



P\ < 



nl/6 



|p(p)r a (i + ipi) ma dp 



(i + ipi)- m6 dp 



1/6 



Thus J |p|* < oo if we can satisfy £a = 2, ma < a and m& > 2, in addition to a -1 = 1. 

This requires a/a > m > 2/6, or 1 < a < 1 + Thus, we require t = 2/a > 4/ (a + 2) 
which, since q = t/(t — 1), means g < 4/(2 — a). Q.E.D. 

Corresponding to the minimization (3.2) there is a variational "equation" satisfied 
by the minimizer p. In the general case, in which p might be a measure, the variational 
equation exists but is slightly complicated to state. 

In physically interesting cases V is certainly in C 1,a , in which case Theorem 3.2 tells 
us that p is a function and that p * is continuous. Hence the total potential 

V c = V + p*\x\~ 1 (3.13) 

is continuous. It is then easy to derive by standard arguments, as in Sect. 1, that p is the 
unique nonnegative solution to the variational equation 

V(x) + p * = \x if p(x) ^ 

V(x) +P* Ix]- 1 > [i if p(x) = 0, (3.14) 

for a unique fx > 0. As usual the chemical potential fi is a monotone function of the particle 
number N = J p. In the special case of a parabolic confining potential the solution to (3.14) 
can be given in closed form. 

3.3 PROPOSITION (Minimizer for the parabolic exterior potential). IfV(x) = K\x\ 2 
then the minimizer of £^ v is 
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where A = (8K/3nN) 2 / 3 . In fact, p% v is the solution to (3.14) with p = (3yr/4)iVA 1/2 . 

Proof: This solution (3.15) was certainly known before; see e.g., [10]. We give the 
proof here for the convenience of the reader. We only have to show that p = P%y is the 
solution to (3.14). It is enough to consider the case A = 1 and N = 1. Then V(x) = ^r\x\ 2 



(\y\ - 


- x| cos6>) 2 


i - 


x 2 sin 2 9 



and \x — We may compute p * |x| 1 = J a/1 — \x — y\ 2 \y\ 1 dy by writing y in polar 
coordinates (|y|,0) and performing the \y\ integration first. 

p*\x\- 1 = ^- J J[(l-\y\-\x\cos9) 2 -\x\ 2 sin 2 6] 1/2 d6d\y\ 

the integrations are over the intervals in 9 and \y\ for which the integrands are real. 
Introducing the variable t = (\y\ — \x\ cos6>)(l — \x\ 2 sin 2 9)~ 1 / 2 we obtain 

p*\x\~ 1 = ^ j (l-\x\ 2 sm 2 9)d9 J {l-t 2 f' 2 dt, 

-0m(x) 



J ri ,.|2,.;,,2 

-0 m (x) 



, , (l-\x\ 2 8m 2 9)d9, (3.16) 



where 



Thus 



and 



Vm[X) - \ gin -l 1 ^ if , > L 



l ) i 37T 3"7T 2 

p * X = X if X < 1 

ii 4 8 1 1 ~ 



_i 37T 37T . l9 
p*\x\ > X if X > 1 



(The last inequality comes from the fact that the integral from 9 m to n/2 in (3.16) is 
negative.) Q.E.D. 
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The energy function E C (N, V) has the simple scaling: 

E C (N,V) = N 2 E C (\, -|Q . (3.17) 
The minimizing density p^ v for (3.2) scales as 

P%, v (x) = Np° v (x). 

where v = V/N. 

We shall now make precise in what sense the classical energy is a limit of the MTF 
energy. In fact, in two different limits (the large B limit and the low coupling limit) the 
MTF energy will converge to the classical energy. We first treat the large B limit. 

3.4 THEOREM (Large B limit of MTF). If the exterior potential V is in the class 
C\oc we have, as B — > 00, 

E MTF (N,B,V) -> E C (N,V) (3.18) 

and 

PnJvW - Pn.v, (3-19) 

in the weak L 1 sense. 

Proof: If we use P^ T B F v as a trial density in £ c and recall that js > Owe immediately 
obtain E C (N, V) < E MTF (N, B, V). 

For the bound in the opposite direction we use P% v as a trial density for £ MTF . In 
order to do this it is, however, important that we know (from Theorem 3.2) that p^ v is a 
function. Hence jB(P%,v) ^ s weu defined. Moreover, from the definition of jb, 3b(Pn,v) ~^ 
almost everywhere as B — > 00 and jB(P%,v) ^ Jo(Pn,v)- Since jo(p) = (7r/2)(p) 2 we 
know from (3.9) that jo(P% v ) is integrable. The limit in (3.18) is therefore an immediate 
consequence of Lebesgue's dominated convergence theorem. 
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The convergence of the densities in (3.19) follows in a standard way by replacing V 
by V + ef with / a bounded (measurable) function and differentiating w.r.t. s, see e.g., 
[29]. Q.E.D. 

We point out that if p^ v is a bounded function, as it is, e.g., for V = K\x\ 2 , then 

E MTF (N,B,V) = E C (N,V), (3.20) 

for B large enough, because in that case Jb{P% y) vanishes for B large. 

Finally, we now discuss the weak coupling limit in the case of homogeneous exterior 
potentials. Suppose V is a homogeneous function of x, V(Xx) = \ s V(x), s > 0. If we 
consider the exterior potentials kV(x) with k > the classical energy and density obey 
the scalings 

E C (N, W) = k 1 ^ s+1 ^E c {N, V) (3.21) 

and 

P%kvi*) = k 2/ ^ P % v (k^x) . (3.22) 

If k is small we see from (3.22) that the minimizing density for the MTF functional 
will spread out and its kinetic energy will be negligible compared with the classical terms. 
We prove this rigorously now. 

3.5 THEOREM (Weak coupling limit of MTF with homogeneous potentials). Let V 
be C^'" and homogeneous of degree s. If k — > then 



E MTF (N,B,kV) 

E G (N, kV) ^ ' 



(3.23) 



and 

fc" 2/(s+1) «V (k-W+»z) - p% v (x) , (3.24) 
in weak L 1 sense. Both limits are uniform in B. 
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Proof: As above we may use P^ T Jy as a trial density in £ c to conclude that E C (N, V) < 
E MTF (N,B,V). 

To prove the bound in the opposite direction we again use P% iV as a trial density for 
^mtf We then ^tain from (2.15) and the scaling (3.22) that 

E MTF (N, B, kV) < k 2 /^ J jo(P%, v ) + E C (N, kV), (3.25) 

where we used that jb < jo- If we compare this with the scaling in (3.21) we see that 
E MTF /E C -> 1 as k -> since jo{p% y ) = (n/2)(p% y ) 2 is integrable. 

The convergence of the densities follows again by replacing V by V + ef and differen- 
tiating w.r.t. e. Q.E.D. 

In the same way as for the large B limit (3.23) becomes an identity for small k if p% v 
is a bounded function. 

We may of course also introduce the scaling V = Nv when v is homogeneous of degree 
s. Then kV = Kv, where K = kN, and the limit in (3.23) is uniform in N. The limit in 
(3.24) is uniform if we formulate it as 

iV-^-^+^V (k-W+Vx) - pl v (x) . (3.26) 

We remark that if a potential W is asymptotically homogeneous in the sense that 
there is a homogeneous potential V with lim^i^oo W(x)/V(x) = 1, then 

lim E(N, B, kW)/k 1/(1+s) = E C (N, V) (3.27) 

k— >0 

uniformly in B, where s is the degree of homogeneity of V. 

IV. THE CLASSICAL POINT CHARGE MODEL: A LIMIT OF QUANTUM 
MECHANICS 



25 



QDOTS EHLJPSJY 3 Fob /94 

Another model that sheds some light on the physics of our problem — and that will 
also be important for bounding the difference between the TF theory and the original 
quantum theory in Sect. V — is the classical particle model. In this model the kinetic 
energy is simply omitted altogether, but the point-like nature of the electrons is retained. 

4-1. DEFINITION (Classical particle energy). With V{x) being the confining poten- 
tial the classical particle energy for N points in R 2 is defined by 

N 

£ p (x u ...,x N ;V) = J2 v ( x i)+ \xi-X3\- 1 . (4.1) 

i=l l<i<j<N 

The minimum classical particle energy for N point particles in R 2 is 

E P (N,V) =m£{S p (x 1 ,...,x N ;V) : x t G R 2 }. (4.2) 

We shall estimate the particle energy E P (N, V) in terms of the classical continuum 
energy E C {N, V). We first show that E C (N, V) gives an exact upper bound on E P (N, V). 
4-2. LEMMA (Upper bound for E p ). For all N we have 

E P (N, V) < E C (N, V) - N 3/2 /(8R v ), (4.3) 

where R v is the maximal radius given in Lemma Al. 

Proof: First, let us give a very simple argument that yields an error term proportional 
to N instead of iV 3 / 2 . The energy E P (N, V) is bounded above by 

J £ P (x 1 , . . . , x N ; V)$(xi, . . . , x N )dx! . ..dx N , 

N 

for any nonnegative function $ with J $ = 1. We take $(xi, . . . ,xn) = Yi ^[l]^)) where 

i=i 

for simplicity we have introduced the notation p^j for the minimizer p^ v i N for £ c [p; V/N] 
with J Pi V / N = 1. Note that p^ depends only onti = V/N . We obtain 

r N 

E P (N,V) < / E p (x 1 ,...,x N ;V) 1 [[p [1] (x i )dx 1 ...dx N 
J i=i 

= N J V(x)p [1] (x)dx+ N ( N 2 ^ J j Pmix^x - yl^Pm^dxdy. 
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C : n C 



Recalling that the minimizer of 8 is P%,v( x ) = ^P[i]( x )i we get an error term — a A, with 
a = h S S P[i]( x )P[i](v) \x -yl^dxdy. 

Now we turn to a proof of (4.3) which, obviously, has to be more complicated than 
the previous discussion. By Lemma Al there is a fixed square, Q, centered at the origin, 
whose width, W, equals 2R V , such that the minimizer p = P%y for S c is supported in Q. 
For simplicity we suppose that y/~N is an integer; if this is not so the following proof can 
be modified in an obvious way. 

First, cut Q into \/~N vertical, disjoint strips, Si, S2, ■ ■ ■ , S ^ such that j s P = VN 
for all j. Let tj denote the width of Sj, so that Y^j^i^j = W. Next, make y/~N — 1 
horizontal cuts in each Sj so that the resulting rectangles Rjk for k = 1, ... , VW satisfy 
f R k p = 1. Denote the height of these rectangles by hjk, so that J2k=i hjk = W for each 
j. Having done this we note, by convexity, that for each j 



;v- 1/2 J> + M" 1 > 



k=i 



[tj + N-Ww]' 1 . 



k=i 

Again, using the same convexity argument for the j-summation, we have that 



N 3/2 

EEfe + M _1 >^ (4-4) 

j=l fc=l 

Let Pjk be the minimizing density p restricted to the rectangle Rjk, i.e. Pjk( x ) = 
1 if x G -Rjfc and = otherwise. Thus, J Pjk = 1. We denote these A function by 

JV 

p\ 2 = 1, . . . , A. Define &(xi, . . . ,xjv) := n and, as in the previous proof, a simple 

computation yields 

r N 

E P (N, V)< J S p ^ = S c [P;V]-J2d(p\P 1 ), 



i=i 



with D(f, g) = UI f( x )9{y)\x - y\ Maxiy. 
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To complete our proof we note that as long as x and y are in Rjk we have that 
\x - y\~ x > (tj + hjk)- 1 . Thus, Y. D (Pjk,Pjk) > N s / 2 /4W by (4.4) and the fact that 
J p jk = 1. Q.E.D. 

4.3. LEMMA (Lower bound for E p ) . Assume that V is a potential in C lo ' c . Then for 
all N we have 

E P (N,V) > E C (N,V) -bN s/2 . (4.5) 

with 

b=isH J (pf,W) V2 dx+ (2R v )- 1+ ' 2 "> (/ «„)') "* (4.6) 

and where q is any number satisfying 2 < q < 4/(2 — a), R v is the maximal radius given 
in Lemma Al and p = q/(q — 1) < 2. As we explained in Theorem 3.2, our hypothesis 
that V E C^'" implies that for q < 4/(2 — a), (J (Pi v ) q ) ^ is less than some constant that 
depends only on q, a and v = V/N . In particular, b depends only on v. 

We see from (4.3) and (4.5) that when V G C lo '" the power 3/2 in the error term is 
optimal. 

In order to prove (4.5) we need the following electrostatics lemma of Lieb and Yau 
[26]. The original version was for R 3 ; we state it here for R 2 solely for the convenience of 
our present application. 

4-4- LEMMA (The interaction of points and densities). Given points x\,...,xn in 
R 2 , we define Voronoi cells Ti, . . . ,Tn C R 2 by 

Tj = {y E R 2 : \y - Xj\ < \y - x k \ for all k ^ j} . 

These Tj have disjoint interiors and their union covers R 2 . We also define Rj to be the 
distance from Xj to the boundary ofTj, i.e., Rj is half the distance of Xj to its nearest 
neighbor. Let p be any (not necessarily positive) function on R 2 . (In general, p can be 
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replaced by a measure, but it is not necessary for us to do so.) Then (with D(f,g) = 
\ II f(x)g(y)\x - y\^dxdy ) 

N f 

E \^~ x j\ 1 >-d(p,p) + J2 / p(y)\y-xj\ 1 d y 

l<i<j<N j=l ^ R2 

N N 

+ i£ i *7 1 -£/ p(y)\y-^\-^y ■ (4-7) 

.7 = 1 3 = 1 JF i 



Proof of (4-5): We choose p in (4.7) to be the minimizer p^ y for the functional 
and we choose the x^s to be any (not necessarily minimizing) configuration for £ p . It is 
important, for us that a minimizer exists for £ c , for then p satisfies the Euler-Lagrange 
equation (3.14). Since J p = N, we conclude from (3.14) that 

N r r i N r r 

E J Piy^y-xjl^dy + Vixj) =Y J V c (x J )>N l i = J V c p = 2D(p,p) + J Vp . 



3 = 1 



J = l 



Thus, if we add J2j V( x j) to both sides of (4.7) we have that 



N 



E p (N,V)>E c (N,V)+J2 l s R J 1 -f P(y)\y ~ xjl-'dy . 

3=1 Jr > 

Our goal will be to control the rightmost term in (4.8) by the RJ 1 term. 

We split each region Tj into two disjoint subregions, 1^ = Aj U B d , where 

Aj := {x : \x — Xj\ < Rj}, Bj := {x G Tj : \x — Xj \ > Rj} 

Then, by Holder's inequality 

N / \ 2/3 / N 1/3 

E/ ly-^-rV(y)dy< W p" /2 \ E/ \v-*i\- z *v 

j=\ jB i \ j J Fj J y j ■) y •<•/ Hi 



(4.8) 



,3/2 



. V3 

2/3 ' 



27rE^7 1 • ( 4 - 9 ) 



3 
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If we define X := V . R- 1 we can rewrite (4.9) plus the X/8 term in (4.8) as — 

2 /3 

(J> 3 / 2 ) 7 X 1 / 3 . The minimum of this quantity, over all values X, is (A/3)y/2/3 J P 3/2 , 
and thus we have accounted for the first error term in (4.6). 



To estimate the term / := XljLi Ia P{v)\v ~ x j\ 1< ^V some control is needed over the 
possible singularities of p. Let p = q/(q — 1) be the dual of q. Then 



r 1 1 / p 



j 



(4.10) 



We note that, since 1 < p < 2, 





i/p 


JV 


l/p-l/2 




< 










J=1 





by Holder's inequality. Now 7ri? 2 is the area of the disc Aj and thus it i? 2 is the total 
area of all these disjoint discs. How large can this area be? To answer this we recall 
Lemma Al in the appendix which states that for the purpose of finding a set of points 
that minimizes the classical particle energy S p we can restrict attention to a disc of radius 
R v , centered at the origin. We may therefore assume that our x/s satisfy \xj\ < R v . This 
we can do whether or not an energy minimizing configuration exists. Having done so and 
assuming that iV > 2 we have that Rj < R v for all j, and hence all our discs are contained 
in a disc of radius 2R V centered at the origin . Thus ^ - R? < (2R V ) 2 , and our second 



error term, (4.10), is bounded above by (/ P q ) /q (^) 1 /p(2R v )- 1+2 /pN 1 / 2 . This yields 
(4.6). Q.E.D. 



V. MTF THEORY IS THE HIGH DENSITY LIMIT OF QUANTUM ME- 
CHANICS 

In this section we prove that the quantum energy and the quantum density are given by 
the corresponding MTF quantities to leading order for large N. These are the statements 
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of Theorems 1.1 and 1.2. We shall not prove Theorem 1.2 since it follows from Theorem 1.1 
in a standard way by replacing V by V + ef with / a function in C lo '" and differentiating 
w.r.t. £, see e.g., [29]. 

We shall prove Theorem 1.1 by giving sharp upper and lower bounds to the quantum 
ground state energy. The upper bound is obtained by a variational calculation using the 
magnetic coherent states introduced in [25,21]. The lower bound is more difficult. Besides 
the results of the previous sections several ingredients are needed. The first ingredient is 
a kinetic energy inequality of the Lieb-Thirring type [31-32], [21], [33]. Such an inequality 
estimates the kinetic energy of a many-body wave function from below in terms of a func- 
tional of the density. The proof of this inequality in the 2-dimensional case considered here 
is harder than the 3-dimensional case treated in [21] and involves some new mathematical 
ideas. The second ingredient is a lower bound on the exchange-correlation energy. The 
proof of this inequality is similar to that given in [34-35] for the 3-dimensional case. 

Once the kinetic energy and exchange-correlation inequalities have been established 
the proof of the lower bound is completed by a coherent states analysis. 

We start by discussing the magnetic coherent states used in the proofs of both the 
upper and lower bounds. They are constructed from the kernels 

U aa (xa\ ya") = ^- exp{z(x x y) ■ B - \x - y\ 2 B/A}L a (\x - y\ 2 B 12)8^8^ (5.1) 

of the projection operators onto the Landau levels a = 0, 1, 2, . . . with ^-component of spin 
a = ±1/2. Here L a are Laguerre polynomials normalized by L a (0) — 1. In fact, all that 
matters are the projectors 11^ on the states with energy e v {B); these are given by a sum 
of at most two of the projections H aa . More precisely, 

a, a 

c C +^+ 1 a=e„(B)/B 
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We shall not need the explicit form (5.1). The three important properties of 11^ that 
we use are the following 

^2 ^(W, ya") =S(x- y)5 a > a », (5.3) 
J2^A^\xa')=d u (B) (5.4) 

a' 

H^Jl v = e v {B)U v , (5.5) 

where ifkin is given by (2.1). 

Let g be a real continuous function on R 2 , with g(x) = for \x\ > 1, f g 2 = 1 
and j{Vg) 2 < oo. (The optimal choice that minimizes j{Vg) 2 is the Bessel function J , 
suitably scaled and normalized.) Define g r (x) = r~ 1 g(x/r), with < r < 1 to be specified 
later. For each ttGR 2 ,z/ = 0,l,2, we define the operator U vu — the coherent "operator" 
— with kernel 

n.uK, ycr") = g r (x - u)U v (xa', ya")g r (y - u). (5.6) 

It easily follows from (5.3,4) and the properties of g that these kernels satisfy the coherent 
operator identities ([36]) 

J IWzct', ya")du = S(x - y)<W" (5.7) 

TrLU = [ ^u{xa',xa')dx = d„{B). (5.8) 
Moreover, a simple computation gives, using (5.5), 

Tr[tf kin rU] = d v (B)[e v (B) + J (Vg r )% (5.9) 
Tr[VU vu ] = d v {B)V * g 2 (u), (5.10) 
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where V is a (continuous) potential and * denotes convolution. Likewise, for all / with 

</!/> = i 

(/|#id»|/> = J fr„(5)</|IU|/)d« - J (Vg r ) 2 (5.11) 

(f\V*9 2 r \f) = J2 \ n«X/|IU|/>d«. (5.12) 

Equations (5.9,10) will be used in proving the upper bound, while (5.11,12) are needed for 
the lower bound. 

5.1 THE UPPER BOUND 

We use the variational principle of [37] . According to this principle 

^N,B,V)<^W^ + V)K\ + \Y. f / ""'"^'^ dri,, (5.13) 

^ Jr 2 Jr 2 \ x -y\ 

for all operators K, with kernel K(xa,ycr'), satisfying 

0<(f\K\f) <(/|/), (5.14) 

for all /, and 

Tr[K]=V / K(xa,xa)dx = N. (5.15) 
^ a Jr 2 

We shall choose K as follows. Let p MTF be the MTF density, i.e., the minimizer of the 
functional (1.5) with J p MTF = N. Denote by z/ max (x) the highest filled level. Then 

0<p MTF (x)- Yl d v {B)<d v _ {x)+1 {B). (5.16) 

We introduce the filling factors 

1 V < ^max(^) 

fAx)= )[p MTF (x)~ £ d„(B)]/d„_ (x)+1 (B), u = u nuac (x) + l (5J?) 
v > z/ max (x) + 1 
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and define 

K(xa,ya') = / fu(u)n vu (xa, ya')du, (5.18) 

with Ii vu as in (5.6). It follows from (5.7) that K satisfies (5.14) and from (5.8), (5.16) 
and (5.17) that Tr[K] = jp MTF (u)du = N. 
Note that (5.4) and (5.6) imply 

K(xa, xa) = p MTF * g 2 r {x). (5.19) 

a 

Hence, the last term in (5.13) is D(p MTF * g 2 , p MTF * g 2 ), where the functional D was 
defined in (1.4). By convexity of D we find that 

D(P MTF * g 2 r , P MTF * g 2 r ) < D(p MTF , p MTF ). 

From (5.9), (5.10) and (5.13) we obtain 



IX 



E<Z(N,B,V) < S MTF (P MTF ) + N J(Vg r ) 2 dx + J [V * g 2 (*) - V(x)]p MTF (x)dz 

< E MTF (N, V, B) + Nr~ 2 [ (Vg{x)) 2 dx + N 2 sup [v * g 2 (x) - u(x)],(5.20) 

J \x\<R 

where R = R v is the finite radius and we have written V = Nv. Since v is in C 1,a 

sup \v * g 2 (x) — v(x) \ < (const. )r. 

\x\<R 

We can choose r = tn such that tn — *• and r^ 2 /N — > as N — > oo. This means that 
r^v should go to zero but still be large compared with the average spacing iV -1 / 2 between 
electrons. The optimal choice is of the order r = (const. )N~ 1 ^ 3 . Thus the error 

[E Q (N, B, V) - E MTF (N, B, V)]N~ 2 

is bounded above by a function e^{v) = c + (v)N~ 1 / 3 (independent of B). This finishes 
the proof of the upper bound. 
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5.2 THEOREM (Lieb-Thirring inequality in 2 dimensions). 

Let Ha = |(*V — A) 2 + S • B. (This is the operator H^ n from (2.1) with 7 = 1.) Let 
W be a locally integrable function and denote by ei(W), €2(W), . . . the negative eigenvalues 
(if any) of the operator H = H A — W defined on L 2 (R 2 ; C 2 ), the space of wave functions 
of a single spin-\ particle. Define |W| + (x) = |[|W"(a;)| + For all < A < 1 we 

then have the estimate 

V \ej(W)\ < A" 1 ^- / \W\ + (x)dx+ |(1- A)" 2 f \W\ 2 + (x)dx. 

Proof: For any self-adjoint operator A we denote by iV a (A) the number of eigenvalues 
of A greater than or equal to a. 

Since replacing W by its positive part \W\ + will only enhance the sum of the negative 
eigenvalues we shall henceforth assume that W is positive, i.e., W = \W\ + . We consider 
the Birman-Schwinger kernel 

K E = W 1/2 (H A + E)- X W 1I<2 . 

According to the Birman-Schwinger principle (see, e.g., [38], p. 89) the number Ne(—H) of 
eigenvalues of H below — E is equal to the number N\{Ke) of eigenvalues of K E greater 
than or equal to 1. We find 

rcc poo 

y>;WI= / N E (-H)dE= / N 1 (K E )dE. 
■ Jo Jo 

In order to estimate Ni(K E ) we decompose the Birman-Schwinger kernel into a part 
K E coming from the lowest Landau level and a part coming from the higher levels. If 
n is the projection onto the lowest Landau band these two parts are defined by 

K% = W 1/2 n (H A + E)~ 1 H W 1 / 2 = E~ 1 W 1 ^ 2 U W 1 ^ 2 (5.21) 
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and 

K> = W X ' 2 (I - U )(H A + E)-\I - U )W 1/2 . 

Since n commutes with Ha we have Ke = K E + K E . 

Now we use Fan's theorem [39], which states that if Hi(X) > ^(X) > ... denote 
the eigenvectors of a self-adjoint compact operator X then /j, n+m+ i(X + Y) < /j, n+ i(X) + 
fj> m+ i(Y) for n, m > 0. From this we have Ni(X + Y) < N X (X) + Ni_ x (Y) which, in our 
case, reads as follows. 

N^Ke) < N X (K%) + N 1 _ X (K>), for < A < 1. 

This inequality permits us to consider the two parts of Ke separately. We first consider 
the contribution from the lowest level: N X (K° E ) = N XE (W 1 / 2 U W 1 / 2 ). We get 

poo poo 

/ N x (K° E )dE = / NxsiW^UoW^dE 
Jo Jo 

poo 

=A _1 / N E (W 1/2 U W 1/2 )dE = \- lr Tr(W l,2 Il W 1/2 ) 
Jo 

=A _1 J U (x,x)W{x)dx = A_1 ^y W(x)dx. 

The second part is straightforward. We first notice that Ha(I — IIo) > B(I — Ho)- Hence 
H A (I-U ) > l(H A + ±B)(I-Il ) > i(zV-A) 2 (7-n ). (Note that (zV-A) 2 commutes 
with n .) Since the operator inequality < X < Y implies X~ x < Y~ l we have that 

K> < W 1/2 (I - n )[|(iV - A) 2 + E]-\I - no)^/ 1 / 2 < W 1/2 [±(iV - A) 2 + E\~ X W X I 2 . 

We conclude that N X - X (K~^) < N^Ke), where 

K E = ((1 - A)- 1 ^) 172 [|(iV - A) 2 + E]- 1 ((1 - X)~ 1 W) 1/2 
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is the Birman-Schwinger kernel for the operator H = |[(iV — A) 2 — 3(1 — A) _1 VF]. The 
Birman-Schwinger principle implies that J °° Ni (KE)dE is the sum of the negative eigenval- 
ues of H. An estimate on this quantity follows from the standard Lieb-Thirring inequality, 
i.e., 

jf N x _ x {K>)dE< (0.24)1 J (^JL_ W ( X ^ dx<\{l-\)- 2 j W{x) 2 dx. 

The constant 0.24 can be found as L 1)2 in [40], Eq. (51). It was improved slightly by 
Blanchard and Stubbe [41], see also [42]. In these references only the case A = was 
considered. It is, however, a simple consequence of the diamagnetic inequality (see [43]) 
that the constant is independent of A. Q.E.D. 
The Lieb-Thirring inequality in Theorem 5.2 implies an estimate on the kinetic energy 

N 



in terms of the one-particle density 



/ , -"kin 
J = l 



P<,(x) = N £ •■• £ / \ip(x, x 2 ...,x N ;a 1 ,..., a N )\ 2 dx 2 . . . dx N . 

<7i=±l/2 CT iV=±l/2j^2(JV-i) 

Here ip is a normalized A particle fermionic wave function. 

5.3 COROLLARY (Kinetic energy inequality in 2 dimensions) . Let and p^ be 
defined as above. Then for all < A < 1 we have 

f0, if P^ < A -1 ^ 

T*>{ . • (5.22) 

I H 1 " A) 2 / [p<,(x) - A" 1 ! ] 2 dx, z/p^ > A" 1 ! 

Proof: The inequality in Theorem 5.2 holds for the operator — W if I7I > 1. 
If I7I < 1, however, one should choose n in the proof of Theorem 5.2 as the projection 
onto the levels v = and v = 1 (not only onto v = 0). Equation (5.21) is then no 
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longer an identity but a bound. In this way one concludes that the negative eigenvalues 
ei(W), e 2 (W), ... for H kin - W satisfy 

y2\ej(W)\<ct [ \W\ + (x)dx + (3 [ \W\ 2 + (x)dx, 
j Jr. 2 Jr 2 

with a = A -1 ^ and (3 = |(1 — A) -2 . This bound is clearly valid for all 7. 

The proof of (5.22) now follows by a standard Legendre transformation. In fact, if 
W > we have 

T^ = (^J2 H kil ~ W ( x j) + J Wp^ > [Wp^ -aW- [3W 2 ] (5.23) 

Since the Legendre transformation of the function W aW + f3W 2 is the function 

P^ sup \pW-aW-l3W 2 ] = I ,_ u , 2 li /- a , 
w>o l J I (4/3) >-«) 2 , ifp>a' 

we see that (5.22) follows by making the optimal choice for W in (5.23). Q.E.D. 

N 

5.4 LEMMA (Exchange inequality in 2 dimensions). Let tp G (^)L 2 (R 2 ;C 2 ) be any 
normalized N particle wave-function (not necessarily fermionic) and let 

N 

Pj>{x) = ^ ^ • • • ^ / \^(x 1 ,...,x N ;a 1 ,...,a N )\ 2 dx 1 ...dx i _ 1 dx i+1 ...dx N 

i=l ai=±l/2 crjv^il/^ajv 

6e £/ie corresponding one-particle density. Then 

E-E /w 2 £ 

— R 2iV l<i<j<N 

>\J J P^(x)p i} (y)\x-y\- 1 dxdy- 192(27r) 1 / 2 J p^(x) 3/2 dx. (5.24) 

R 2 R 2 R 2 



Proof: The proof is essentially the same as in [34] , where the three dimensional equiv- 
alent of (1.1) was proved. Our presentation is inspired by [44]. 
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We use the representation (in 3 dimensions a similar representation was originally 
used by Fefferman and de La Llave [45]) 



\x-y\- L =n- L J J Xr(x- z)x R (y- z)R-*dRdz, (5.25) 

R 2 R 

where Xr is the characteristic function of the ball of radius R centered at the origin. If we 
use (5.25) to represent Yl \ x i ~ x j\~ X we can estimate the integrand as follows. 

i<j 

X x r( x * - z )Xr{ x j ~ z ) = \ I ^Xnixi - z) J - \^2,XR{xi - z) 

l<i<j<N V i J i 

Pip{y)X R {y- z)dy^ -\^X R {xi-z) 
>^X R {xi-z) j P^(y)X R (y- z)dy-\(^J ' P^{y)X R {y - z)dy^j -\^X R {xi- z). 
If we integrate this inequality over the measure R~ A dRdz, the last term, | Yl XR(xi — z), 

i 

will give a divergent integral. For the purpose of a lower bound, however, we can restrict 
the integration in (5.25) to R > r(z), where r(z) > is some specific function we shall 
choose below. Using the fact that 

53 •••X) / \^\ 2 J2x R (x l - z)dx 1 ...dx N = J p^(y)x R (y - z)dy, 

we obtain, 

^ • • -XI J H 2 X \xj - xjl^dx! . ..dx N 

J Pip{y)x R {y - z)dy^j R~ 4 dRdz 



> K 1 




R>r(z) 
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1^-1 

2 71 



R 2 

> 




R>r(z) 



Pi>(y)X R (y - z)dyR A dRdz 



\tt 1 J J P^{x)p^{y)\x - y\ 1 dxdy 



-iyr- 1 / / 7r 2 p^(z) 2 dRdz- / / R- 2 Kp*^{z)dRdz. (5.26) 

J JR<r(z) J J R>r(z) 

R 2 R 2 

Here we have introduced the Hardy-Littlewood maximal function 

P^(z) = sMnR 2 )- 1 J P^(y)x R (y - z)dy, 

which, viewed as a map from L P (R 2 ) to L P (R 2 ), is a bounded map for all p > 1 (see [46], 
pp. 54-58). The error terms in (5.26) can be computed as 



K 1 



R 2 



7i 2 P* i ,(zYdR+ / R- 2 7tp^(z)dR 

R<r(z) JR>r(z) 



dz 



= K 1 / (r(z)n 2 p;(z) 2 + r{z)-\p^{z))dz. 



R 2 

The optimal choice for r(z) is r(z) = (np*(z))~ 1 ^ 2 . This means that the error is 7T 1 / 2 J P^(z) 3 / 

R 2 

but this can be estimated by the maximal inequality to be less than 192 (27T) 1 / 2 J p(z) 3 / 2 dz. 

R 2 

Q.E.D. 



5.5 THE LOWER BOUND 

Our goal here is to give a lower bound to i? Q (iV, B, V) in terms of E MTF (N, B, V) with 
errors of lower order than N 2 as N tends to infinity. It is important here that V = Nv, 
where v is fixed. To be more precise we shall prove that 

N- 2 [E^(N, B, V) - E MTF (N, B, V)] > -e^(v), (5.27) 

where e^(v) is a non-negative function which tends to as N — > oo for fixed v. Note, 
however, that sj^(v) does not depend on B. 
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We shall treat the cases of large B and small B separately. In the large B regime we 
prove (5.27) by a comparison with the classical models discussed in Sects. Ill and IV. In 
the small B regime we use magnetic coherent states, Theorem 5.2 and Lemma 5.4. 

The dividing line between large and small B is determined as follows. If the minimizer 
Pi v of S c , with J Pi v = 1 and confining potential v, is bounded (e.g. for v(x) = \x\ 2 ) 
then we define small B to mean 

B/N < p c := 2n sup p^ v (x) (5.28) 

X 

As explained in (3.20) we have for (3 > (3 C that £ MTF (1, 0, v) = E c (l,v). 
For the general class of v where we do not know that the minimizer is bounded 
we simply define 

Pc = N 1 / 3 . 

By Theorem 3.4 we then have that the function 

S(N,v)= sup \E MTF (l,[3,v)-E c (l,v)\ (5.29) 

tends to zero as iV tends to infinity. 

Case 1, B/N > (3 C : By simply ignoring the kinetic energy operator, which we had 
normalized to be positive, we have the obvious inequality E Ci (N, B, V) > E P (N, V) where 
E p is the energy of the classical point problem. 

From Lemma 4.3 we can therefore conclude that 

E Q (N, B, V) > E P (N, V) > E C (N, V) - b(v)N 3/2 . (5.30) 

Since E C (N,V) = N 2 E c (l 7 v) and E MTF (N,B,V) = N 2 E MTF (l, B/N, v) we have from 
(5.29) that 
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E Q (N, B, V) > E MTF (N, B, V) - 8(N, v)N 2 - b(v)N 3/2 . (5.31) 

Thus (5.27) holds with e~ (v) = S(N, v) + b^N' 1 / 2 . 

We emphasize again that if is bounded (e.g. for v = k\x\ 2 ) then 5(N,v) is not 
needed. 

Case 2, B/N < j3 c : In this case we use inequality (5.24) to reduce the many-body 
problem to a one-body problem. 

Let ip be the many-body ground state* for H^. The correlation estimate (5.24) gives 

N 

E^(N,B,V) = (i;\H N \^)>Y,(m^ + V(^m + D(p i} ,p i} )-C / p S J\x)dx. (5.32) 

7 = 1 Jo 



R 2 



We first estimate the last term in (5.32) in terms of the kinetic energy = ( ip 
of ip. According to (5.22) we have 



N 

V H {j) 

2^ -"kin 
J = l 



1/2 / r \ V 2 

2 



J P 3 J 2 < (const.)S 1 / 2 J Pqp + (const.) Q P^j ( f P 



P^,<(const.)B /9^> (const. )B 



< (const. )(t/WcN 3/2 + y%N l/2 ). (5.33) 
Hence, for all < e < 1 (we shall later choose e ~ iV -1 / 2 ) we have 



£Q(iV, B, V) > YfWtt ~ e)H^l + V( Xj )\iP) + D(p^ Pi> ) - (const.) {^W C N 3 / 2 + e^N). 

(5.34) 



where we have used that 



eT^ - (const.) ^/T^N 1/2 > -(const.)£ _1 iV. (5.35) 



Since the exterior potential V tends to infinity at infinity Hn will have a ground state. 
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To relate (5.34) to the MTF problem we use the inequality 
0< J D(^-P MTF ,^-P MTF ) = - 



(5.36) 

which is a consequence of the positive definiteness of the kernel \x — y\~ x . Inserting this 
in (5.34) gives 



N 



E*{N, B, V) > J>|(1 - e)H£l + V(x 3 ) + p MTF * \ Xj \~ V) 

3 = 1 

- D(p MTF , P MTF ) - (const.)( v^iV 3/2 + e^N). 



(5.37) 



Since we have normalized the potential to be positive we have that (1 — e) 1 V(x) > V(x) 
and also (1 - £ )"V MTF * M" 1 > P MTF * M" 1 - 
Hence 

E*{N,B,V)>{l-e) ^E^in + ^O+P™*!^!- 1 )^ 

- D{p MTF ,p MTF ) - (const.) (v^iV 3 / 2 + e - 1 N). (5.38) 



TV 



Obviously 

V { H ^l + V M) +P MTF * N" 1 A > £e„ (5-39) 

3=1 I j=l 

where ei, €2, ■ ■ ■ , e^v are the iV lowest eigenvalues of the one-particle Hamiltonian 



MTF 



^ kin + V(x) + P MTF * Ix]- 1 = H kin + V M '" F (x) 



.MTF, 



(5.40) 



We shall estimate j=i e j by a straightforward coherent states analysis. 

Let /i,...,/iv be the iV lowest normalized eigenfunctions of i/J vITF . For technical 
reasons we introduce a modified operator ff] VITF which is obtained from i7| vITF by replacing 
yMTF ky truncated potential 

T>MTF W = /V MTF W \X\<R V 

K > \CN \x\ > R v ' 
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where R v is the finite radius given in the appendix and C = inf\ x \ >Rv V MTF (x)/N is 
independent of N by the scaling (2.9) of MTF theory. Note that V MTF > V MTF . Then 
from (5.11,12) we have 



N ,N N 



r ~ N 

= J2j du(e v (B) + V*™( U ))'£(f j \IL l/n \f j 



3 = 1 

N 



-r~*N [(Vgf + J2(W MTF - V M ™ tgllfj). (5.41) 

J 3 = 1 

We first consider the last term. Writing \fj(x)\ 2 = p(x) we have 

3=1 

Yf{f 3 \V^ F - V™™ * g 2 r (x)\fj) = J [V MTF (x) - F MTF * g 2 (x)]p(x)dx 

J = 1 \x\<R v +r 

> J [V MTF (x) - V MTF * g 2 (x)]p(x)dx - J V MTF * g 2 r (x)p(x)dx. 

\x\<R v —r R v — r< \x\ <R v +r 

Since V MTF = Nv MTF and jp = Nwe have 



(/,-|F MTF -y MTF *^|/.) > -N 2 sup \v MTF (x)-v MTF *g 2 (x)\-( sup ^ MTF (x))i? 

j =1 \x\<R v \x\<R v +r 

(5.42) 

We can then write (5.41) as 

N r n 

E e ^E/M S )+ yMTF ( w ))E (f3\K„u\fi)du-N 2 e N {v), (5.43) 

3=1 v 3=1 

where 



e N (v)= sup | V MTF (:r)-<; MTF *<7 2 0r) 

|x|<.R„ 



+ ( sup v MTF (x))R 2 v r _ at-V 2 / (Vs) 2 

< C^iV" 1 / 3 . (5.44) 
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For the last step we made the choice r ~ A^ -1 / 3 . 

We focus next on the first term in (5.43). It has the form 



(e u (B) + V MTF (u))p u (u)du : (5.45) 



JV 

where we have denoted Yl ifj\^u\fj) by P v {u). These functions satisfy 



3 = 1 



< p v {u) < TrU vu = d v {B) (5.46) 



and 



J2 u J Pu{u)du = N. (5.47) 

We obtain a lower bound to (5.45) by minimizing over all functions p v satisfying (5.46) 
and (5.47). 

Minimizers p v can be constructed as follows. There is a p > such that 

(d v (B) ife„(B) + V MTF (u)<p 
Pv{u) = Iq ife u (B) + V MTF (u) >p (5.48) 

[<d v (B) \ie u {B) + V MTF {u) = ^ 

All families p v satisfying (5.48) and the constraint (5.47) are minimizers. Note that it 
is possible that e u (B) + V MTF (u) = \x on an open set of u values. The minimizers are 
therefore not necessarily unique. The chemical potential p is uniquely determined by (5.48) 
and the condition (5.47). 

We shall now prove that p = p MTF . All we have to show is that we can find functions 
P u satisfying (5.47) and (5.48) with p = p MTF . 

We know from the MTF equation, Theorem 2.2, that if p MTF (u) = then V MTF (u) > 
p MTF . Since V MTF differs from V MTF only on the set P MTF = we may in (5.48) when 
H = /j MTF replace V MTF by V MTF . We then know from the MTF equation that there are 
unique functions p u satisfying (5.48) with p = p MTF and P u (u) = p MTF (u). In fact, in 

V 
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terms of the filling factors (5.17) we have p v {u) = f v {u)d v {B). Since / p MTF (u)du = N 
we have produced the functions p u allowing us to conclude that p is indeed equal to // MTF . 
If we insert these functions in (5.45) we obtain 

£ J e„(B)p„(«)d« + J y MTF (n)p MTF (n)d W = |(j B (p MTF ( W )) + y MTF ( W )p MTF ( W ))d W , 

(5.49) 

where the identity follows from (2.3). 

We can now combine (5.38), (5.39), (5.43), (5.44) and (5.49) to arrive at 

£ Q (iV, B, V)>(1- e)E MTF (N, B, V) - eD(p MTF , p MTF ) 

- (const.) (^/P~ C N S/2 + £ _1 iV) - C(v)N 5 / 3 (5.50) 

Hence, since D(p MTF , p MTF ) < E MTF (N,B,V) we have 

N~ 2 (E ( ^(N 1 B, V) — E MTF (N, B, V)) > -2eE MTF (l, B/N, v) 

- (const.) (y^WcN~ 1/2 + e^N' 1 + c(v)N~ 1/3 ). 

Note that E MTF (l, B/N,v) is bounded by a constant depending only on v. If we choose 
e ~ iV -1 / 2 we find 

N- 2 (E Q (N,B,V) -E MTF (N,B,V) > -c~ (vfty/fcN-V 2 + N' 1 / 3 }. 

This is equivalent to (5.27) with e^(v) = c~ (f)[ v / ^A r " 1/2 + iV~ 1/3 ]. In the case when 
p c is bounded (3 C is a constant, otherwise we chose it to be (3 C = iV 1 / 3 . In both cases will 
ejf(v) tend to zero as N tends to infinity. This finishes the proof of the lower bound. 

We have proved (1.7). The limits in (1.8) and (1.9) follow immediately from the 
corresponding results for E MTF proved in Sects. II and III. 

5.6 HOMOGENEOUS EXTERIOR POTENTIALS 
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Finally, we shall show how to prove the stronger result Theorem 1.3 for homogeneous 
exterior potentials. 

In this case we do not have that the minimizing MTF density is supported within a 
fixed ball. In fact, the density will spread out as the coupling constant becomes small. 

We shall prove that given e > and ko there is an N £ independent of B such that for 
N > N e and K/N < k 

\E Q (N,B,Kv)/E MTF (N,B,Kv)-l\ < e. (5.51) 

We consider large and small K in very much the same way as we did for B in the 
lower bound above. We shall see below that we can find a k c (depending only on v and e) 
such that (5.51) holds for K/N < k c . 

In the case of large K, i.e., K/N > k c (but K/N < k ) the proof of (5.51) is then 
identical to the proof of Theorem 1.1 given above. 

For small K we again consider the upper and lower bounds separately. We begin with 
the upper bound. We proceed as in Sect. 5.1. We define the trial operator as in (5.17) and 
(5.18) except that we replace p MTF by p% 

Kv The estimate (5.20) now becomes 
E^(N,B,Kv) < £ MTF [p% Kv ;B,Kv]+Nr- 2 [ (X7g(x)) 2 dx + NK sup [v*g 2 {x) -v{x)}. 

J \x\<R K 

Here Rk is the radius of the ball containing the support of p% KV for K = kN. According 
to (3.22) Rk = /c _1 /( s+1 ).Ri, where R\ is the radius for k = 1, which depends only on v. 
Using the homogeneity of v we have 

v * g 2 (x) - v(x) = \x\ s / [v ((x - y)/\x\) - v{x/\x\)\ g 2 (y)dy < c 3 (w)r|x| s " 1 . 
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Thus from (5c. 1) we obtain since s > 1 



E q (N, B, Kv) - S MTF [p c ; B, Kv] < c±(v)N [r~ 2 + KrR s ~ l ] 



< c 5 (v)N r~ 2 + Kr(K/N)-^-^^ s+1 ^ 



<c 5 (v)N 2 / 3 (K/N) 4 /W s+1 », 



with the choice r = (if/iV)~ 2/(3(s+1)) . We also know from (3.25) that 




Hence, from (3.21) and (3.17) we obtain 



E Q (N,B,Kv) < 1 + 



< 




iV- 1 /3(K/ iV )4/(3( s +i)) + (tf/jv) 2 /(*+i)j J E c (N,Kv) 
K/Ar)V(3( s +i)) + ^K/Ny^+^Y) E C (N, Kv) 



(recall that E (N, Kv) is negative, hence the minus sign in the second line). It therefore 
follows from Theorem 3.5 that we can find k c depending only on e (but not on B) such 
that E^(N, B, Kv) /E MTF (N, B, Kv) > 1 - e for K/N < k c . 

We turn next to the lower bound. As in Sect. 5.5 (in the case B/N > f3 c ) we may 
ignore the kinetic energy operator, which we had normalized to be positive. We then have 
the obvious inequality E Cl (N, £?, V) > E P (N,V) where E p is the energy of the classical 
point problem. 

We shall use Lemma 4.3 to compare E p (N,Kv) to E C (N, Kv). We must, however, 
first discuss the scaling of E F (N, Kv). It is clear that if v is homogeneous of degree s then 



£ p (x 1 ,...,x 1 ; Kv) = k 1 /^£ p (k 1 /^x 1 , . . . , k^+^x^v). 
Therefore, E P (N, kV) = ^/(s+i) e p (N, V), i.e., E P (N, kV) has the same scaling as E C (N, Kv) 



[see (3.21)]. 
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From Lemma 4.3 we thus find that 

E p (N,Kv) > E c (N 1 Kv)-b(v)(K/N) 1 ^ s+1 ^N s/2 > E C (N, Kv)(l+cs(v)N~ 1 ^ 2 ). (5.52) 

According to Theorem 3.5 we may thus assume that k c is such that 

E Q (N, B, Kv)/E MTF (N, B, Kv) < (1 + e/2)(l + c s {v)N- 1 ' 2 ) (5.53) 

for K/N < k c . We can therefore clearly find N e such that the right side of (5.53) is less 
than 1 + e for N > N e . 

VI. CONCLUSION 

We have analyzed the ground state of a two-dimensional gas of iV electrons interacting 
with each other via the (three-dimensional) Coulomb potential and subject to a confining 
exterior potential V(x) = Kv(x) where K is an adjustable coupling constant. The electrons 
are also subject to a uniform magnetic field B perpendicular to the two dimensional plane. 

We have found the exact energy and electron density function p(x) to leading order in 
I/AT, i.e., in the high density limit. This limit is achieved by letting K be proportional to N 
as N — > oo, thus effectively confining the electrons to a fixed region of space, independent 
of N. 

It turns out that the answer to the problem depends critically on the behavior of B 
as iV — > oo. There are three regimes. 

(i) If B/N — > 0, i.e., iV 3> B in appropriate units, then normal (two-dimensional) 
Thomas-Fermi theory gives the exact description. Correlations can be ignored to 
leading order in this high density situation. 

(ii) If B/N =constant, a modified TF theory in which the "kinetic energy density" is 
changed from (const. )p 2 to a certain Independent function of p (called Jb{p)) is exact. 
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(iii) If B/N — > oo then the kinetic energy term can be omitted entirely and a classical con- 
tinuum electrostatics theory emerges as the exact theory. This electrostatics problem 
is mathematically interesting in its own right and can be solved in closed form for the 
customary choice v(x) = \x\ 2 . 

Related to the continuum problem is an electrostatics problem for point charged parti- 
cles. Apart from its mathematical interest, it provides a crucial lower bound to the energy 
in case (iii). Another technical point of some interest is the extension of the Lieb-Thirring 
inequality to the two-dimensional particles in a magnetic field which involves dealing with 
a continuum of zero energy modes (i.e., the lowest Landau level). 

APPENDIX 

Here we prove that the minimizers for our three semi-classical problems can be sought 
among densities that vanish outside some finite radius - - for which we give an upper 
bound. This lemma is in an appendix because it pertains to several sections of the paper. 

A.l. LEMMA (Finite radius of minimizers). Consider the three cases: (a) The 
classical energy; (b) the classical particle energy; (c) the MTF energy. Let V(x) be the 
confining potential. We assume that V(x) — > +00 as \x\ — ■> 00 in the sense that the number 
W(R) := inf{V(a;) : |x| > R} tends to 00 as R — > 00. 

Then there is a radius R v , depending only on v = V/N 

E P (N, V) = inf {£ p (x u ...,x N ):\ Xi \<R v for all i) 

E C (N, V) = inf j£ C [/x; V] : support p C {x : \x\ < R v }, J dfj, = AT j 

£ MTF (iV, B, V) = inf |^ MTF [P; B, V] : p(x) = for |x| > R v , J P = Ar| . (a.l) 

Furthermore, any minimizing particle distribution measure or density satisfies the condi- 
tions given in braces in (a.l) 
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A choice for R v , which is far from optimal, is any R satisfying the inequality 

±W(R)>(2 + n- 1 ) + ±(V} 1 , (a.2) 
with (V)i being the average of V in the unit disc: 

(V)i = - f V{x)dx. 

* J\x\<l 

Proof: Particle case: Suppose that \x\\ > R v . Then we move particle 1 inside D, the 
unit disc centered at the origin. The point, y, to which we move particle 1 is not known, 
so we average the energy over all choices of y E D. If we show that this average energy is 
less than the original energy then we know that there is some point y G D such that the 
energy is lowered. Thus, we have to show that 

N ^ 1 N 

Vfa) + J2\xi -xj\ 1 >("^)i + -V/ \y~Xj\ x dy. 

3=2 71 3=2 JD 

Noting that f D \y — x\~ l dy < f D \y\~ 1 dy = 2-7T, by a simple rearrangement inequality, we 
see that it suffices to have W(R V ) > (V)i + 2N, which agrees with (a.2). 

The Classical case: If fx is any measure with J d/i = N, we define (i + to be fx restricted 
to the complement of the closed disc of radius R v centered at the origin. Thus fi + (A) = 
li(A fl {x : \x\ > R v }). Similarly, yT is fi restricted to the disc, so that fi = (i + + yT . 
Assuming that 

fx + 7^ 0, we replace fx by fx e := (1 — e)fx + + fx~ + Sv, where v is Lebesgue measure 
restricted to the unit disc, D, and where rcS = e J dfx + . Thus J dfx e = N. The change in 
energy, to 0(e) as e I 0, is easily seen to be 



5 / V(x)dx — e \ V (x) fx + (dx) + 5 / / \x — y\ 1 dxfx(dy) — e / / \x — y\ 1 fi 
Jd J JdJr 2 J J 



(dx)fx(dy) 
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< 7r<y(V)i - eW(R v ) J dp + + 2n5N, which is negative by (a.l). 

The MTF case: This is similar to the classical case, but with two differences: (i) The 
measure p is replaced by a function p with J p{x)dx = N and (ii) a "kinetic energy" term 
J 3B{P{x))dx is added to the energy. Point (i) only simplifies matters. For point (ii) we 
note the simple fact that jb{p) is bounded above by irp 2 /2 and its derivative, f(p) is 
bounded above by nP; this is true for all B. Let us assume that dp + := p + (x)dx is not 
zero, with P + (x) = p(x) for \x\ > R v and P + (x) = otherwise. The argument is as before, 
but now we must take into account the change in kinetic energy which, to leading order in 
£, is 

$ / JB(P( x ))dx-e / (p + (x))dx < 5n / p(x)dx < SN = e / p + (x)dx. 
Jd Jb Jd J 

The total energy change is then negative by (a.l). Q.E.D. 
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